Astronomy & Astrophysics manuscript no. 3302main 


February 2, 2008 


(DOI: will be inserted by hand later) 





Extinction properties of the X-ray bright/optically faint afterglow 

of GRB 020405 

G. Stratta 1 , R. Perna 2 - 4 , D. Lazzati 2 , F. Fiore 3 , L.A. Antonelli 3 , M.L. Conciatore 3 

1 LATT, Laboratoire d'Astrophysique Toulouse-Tarbes, 14, av. Edouard Belin, 31400 Toulouse, France 

2 JILA and Department of Astrophysical and Planetary Sciences, University of Colorado at Boulder, 440 UCB, Boulder, CO, 
80309, USA 

3 INAF, Osservatorio Astronomico di Roma, via Frascati 33, Monteporzio Catone, 00040, Rome, Italy 

4 Department of Astrophysical Sciences, 4 Ivy Lane, Princeton, NJ, 08544, USA 

Received ; accepted 

Abstract. We present an optical-to-X-ray spectral analysis of the afterglow of GRB 020405. The optical spectral energy dis- 
tribution not corrected for the extragalactic extinction is significantly below the X-ray extrapolation of the single powerlaw 
spectral model suggested by multiwavelength studies. We investigate whether considerable extinction could explain the ob- 
served spectral "mismatch" by testing several types of extinction curves. For the first time we test extinction curves computed 
with time-dependent numerical simulations of dust grains destruction by the burst radiation. We find that an extinction law 
weakly dependent on wavelength can reconcile the unabsorbed optical and X-ray data with the expected synchrotron spectrum. 
A gray extinction law can be provided by a dust grain size distribution biased toward large grains. 
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1. Introduction 

Long-duration Gamma Ray Bursts (GRBs) are associated with 
the core collapse of massive stars exploding as type Ic su- 
pernovae (GRB 980425, Galama et al. 1998; GRB 011211, 
Delia Valle et al. 2003; GRB 030329, Stanek et al. 2003; 
GRB 031203, Malesani et al. 2004). Due to their short life- 
time, massive stars are likely to die within their birthplace. 
Consequently, long-duration GRBs occur in the same star for- 
mation region where their progenitor star was born and where 
it rapidly evolved. If high redshift galaxy star formation regions 
are akin to Galactic giant molecular clouds, a dense and dusty 
environment is expected in the vicinity of GRBs. 

Indeed, the measured equivalent hydrogen column densi- 
ties Nh obtained from X-ray afterglow spectral analysis are on 
average consistent with those observed along the line of sight 
of Galactic molecular clouds (e.g. Galama & Wijers 2001, De 
Pasquale et al. 2003, Stratta et al. 2004). In some afterglows 
(e.g. GRB 971214, Ramaprakash et al. 1998; GRB 980703 
Bloom et al. 1998; GRB 980329, Yost et al. 2002) a high vi- 
sual extinction Ay in the GRB host galaxy has been measured. 
In addition, a large dust content was inferred from high res- 
olution spectroscopy by refractory metal abundance analysis 
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along the line of sight of four optical afterglows (Savaglio et 
al. 2003, 2004). However, on average, the estimated afterglow 
reddening is low and the rest frame visual extinction is a factor 
10-100 lower than expected from the extrapolation of the X- 
ray Nh, assuming the Galactic dust-to-gas ratio (Simon et al. 
2001 ; Galama & Wijers 2000). In a few high redshift GRBs for 
which high resolution afterglow spectroscopy was performed, 
high Nh from Ly a absorption was measured. On the other hand 
a low Ay is inferred from continuum absorption modeled with 
the Small Magellanic Cloud (SMC) extinction law. These re- 
sults can be explained by a low metallicity and/or a low dust- 
to-gas ratio in the burst environment (e.g. Hjorth et al. 2003, 
Vreeswijk et al. 2004). 

Alternatively, the GRB environment may be characterized 
by a "non-standard" extinction law (neither Galactic nor SMC; 
Stratta et al. 2004, Savaglio et al. 2004). In particular, an extinc- 
tion law weakly dependent on wavelength can provide high vi- 
sual extinction without substantial reddening (e.g. Hjorth et al. 
2003). Dust properties of high redshift environments - such as 
the GRB host galaxies - are still poorly known. We have some 
indications from the Magellanic Clouds (e.g. Pei et al. 1992) 
and from a sample of local starburst galaxies (e.g. Calzetti et 
al. 1997) and Active Galactic Nuclei (e.g. Maiolino, Marconi 
& Oliva 2001, Maiolino 2004) that the dust composition, grain 
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size distribution and dust to gas ratio can be significantly dif- 
ferent from the ones observed in our Galaxy. The GRB star- 
forming environment, in addition, may not be representative of 
the typical low-density host galaxy inter-stellar matter (ISM). 
Finally, the intense X-ray and UV flux from the burst itself can 
modify the intrinsic dust grain size distribution by differential 
dust grain destruction, since sublimation processes are more ef- 
fective on small grains (e.g. Waxman & Draine 2000; Fruchter 
et al. 2001; Draine & Hao 2002; Perna & Lazzati 2002; Perna, 
Lazzati & Fiore 2003). 

A diagnostic tool of the GRB environment extinction prop- 
erties is the simultaneous optical-to-X-ray afterglow continuum 
spectral analysis. The intrinsic optical continuum can be ex- 
trapolated from the X-ray data and dust absorption can be mea- 
sured from the deviations of the measurements from the ex- 
pected fluxes. Good quality data are necessary to constrain the 
free parameters of the continuum model (the position of the 
cooling break possibly laying between the optical and X-ray 
bands). This method was applied, with a various degree of ac- 
curacy, in a handful of GRBs (Stratta et al. 2004). 

In this work we present a study of the optically faint/X- 
ray bright afterglow of GRB 020405 (e.g. Dado et al. 2002, 
Masetti et al. 2003, Bersier et al. 2003, Co vino et al. 2003, 
Price et al. 2003, Berger et al. 2003, Mirabal et al. 2003). 
This burst was discovered by the IPN on 2002 April 5.029 UT. 
The duration of the burst was ~ 40s, the 25-100 keV fluence 
was ~ 3 X 10~ 5 ergs cirr 2 and the peak flux was ~ 10~ 6 ergs 
cirr 2 s -1 . This GRB was also observed by the GRB Monitor 
on board BeppoSAX with a duration of ~ 60s in the 40-700 
keV band and a 50-700 keV fluence of 4 x 10~ 5 ergs cirr 2 s _1 
(e.g. Price et al. 2003). Optical observations started 18 hours 
after the burst event (Price et al. 2003) for a period of ~ 10 
days. An unknown fading source was localized at R.A.I 3*58'" 
03U2 and Dec. -31° 22' 22".2 with an uncertainty of 0".3 
(Masetti et al. 2002, Price et al. 2003). An exhaustive sum- 
mary of all the afterglow observations performed in the opti- 
cal band is presented in Masetti et al. (2003). The host galaxy 
has redshift z = 0.691 + 0.002 (Masetti et al. 2003). Chandra 
observations started on April 6.711 UT and lasted until 7.350 
UT using LETGS in conjunction with the ACIS detector and 
revealed a new fading source that was identified as the X-ray 
afterglow counterpart. Mirabal et al. (2003) found a featureless 
spectrum, well described by a power law continuum with en- 
ergy spectral index ax = 0.72 + 0.21 and a rest frame Nh of 
(4.7+3.7)xl0 21 cirr 2 . The X-ray light curve decayed as a power 
law, with a temporal decay index 6x = 1.97 + 1.10 (Mirabal et 
al. 2003). The optical (VRI) decay index is 5 = 1.54 ± 0.06 
from 1 to 10 days after the burst (Masetti et al. 2003). The NIR 
bands (J and H) show a less steep decay with index 6jr ~ 1.3 
(Masetti et al. 2003). Radio observations show a rapidly fad- 
ing "radio flare" at 1.2 days from the burst event (Berger et al. 
2003). Berger et al. (2003) found that the radio, optical and X- 
ray data are self-consistently modeled assuming a collimated 
ejecta expanding into a uniform medium. Best fit parameters 
yield a jet break time at about 1 day after the burst and indi- 
cate that the cooling frequency at that time is at energies lower 
than the optical band. The same conclusion on the location of 
the cooling frequency was obtained by Masetti et al. (2003), 



who showed how the optical spectral energy distribution has a 
break around the J band. For these reasons, for this burst there 
is convincing evidence for the absence of a spectral break be- 
tween the X-ray and the optical band. The different decay index 
and spectral slope observed in the NIR band with respect to the 
optical may be due to the presence of the cooling frequency 
between the optical and the NIR bands (Masetti et al. 2003; 
Berger et al. 2003). 

The paper is organized as follows: in §2 we summarize the 
data analysis procedure; in §3 we describe the adopted extinc- 
tion curves; in §4 and §5 we present and discuss our results. 

2. Data analysis procedure 

We followed the standard procedure for the X-rays (0.1-10.0 
keV) data reduction for Chandra ACIS data. We grouped en- 
ergy channels in order to have at least 20 counts per bin in 
order to apply the x 2 statistic in the fitting procedure. We first 
fit the X-rays spectrum with a power law model and two ab- 
sorber components. The column density of the first absorber 
was fixed to the Galactic value toward the line of sight of the 
burst 4.3 x 10 20 ctrT 2 (Dickey & Lockmann 1995). We subse- 
quently allowed for the second absorber a free column den- 
sity to constrain a possible contribution of absorption from the 
host galaxy ISM. We found that the data are well fitted by this 
model, with a best fit energy spectral index of ax = 1 -0 + 0.2 
and rest frame (z = 0.691) hydrogen column density of Nh z = 
(0.8 + 0.2) x 10 22 ctrT 2 in addition to the Galactic absorption 
0f 2 /d.o.f.=22.3/24, where d.o.f. stands for degrees of freedom). 
The best fit model (absorbed power law) 1.6-10.0 keV flux is 
8.7 x 10~ 13 ergs cirr 2 s _1 . Errors are at 1 cr level. We note that 
our spectral index is steeper than, but still consistent with, the 
one obtained by Mirabal et al. (2003). 

The optical-NIR SED was obtained by taking the photo- 
metric points from the journal table published by Masetti et al. 
(2003). Observations in each photometric band (U,B,V,R,I,J,H 
and K) were performed at different times and with different 
telescopes. We extrapolated the magnitudes at the observa- 
tional epoch to of 7.01 UT (1.98 days after the burst). To 
this purpose, we firstly selected the observations performed as 
closer as possible to f () , and then we extrapolated the magnitude 
at to assuming a single powerlaw with decay index estimated at 
that time from Masetti et al. (2003). The decay index uncer- 
tainty was propagated to the magnitude errors with an average 
increase that depends on the temporal distance of the photo- 
metric measure from to (see Tab.l in Masetti et al. 2003). This 
results in an increase of the error on the magnitudes of 20% 
or more, depending on the temporal distance from the selected 
to- We choose to as the average time of the X-ray observations 
given the large uncertainty in the X-ray decay slope (Mirabal 
et al. 2003). 

We corrected the magnitudes for Galactic extinction (E(B- 
V)=0.054 mag) using the Galactic dust infrared map by 
Schlegel et al. (1998). Assuming a total-to-selective extinction 
Rv = A V /E(B - V) = 3.1, we derived the extinction value 
in each photometric band from the Galactic extinction curve 
parameterization A(A) /A v by Cardelli et al. (1989). We finally 
converted magnitudes to fluxes using the effective wavelengths 
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Table 1. The optical-near infrared magnitudes of the optical afterglow of GRB 020405 (from Masetti et al. 2003), corrected for 
Galactic extinction and rescaled at fo=UT 7.01 (Tgrb+1.98 days) 



U B V R 1 J H K 

21.22 + 0.12 21.90 ±0.1 21.50 ±0.08 21.08 ±0.06 20.56 ±0.12 19.37 ±0.10 18.72 ±0.10 17.99 ±0.10 
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Fig. 1. Extinction curves adopted in the analysis of the optical- 
NIR photometry against 1/A where A is the rest frame wave- 
length. We note that the weakest wavelength dependence is 
shown by the Q2 extinction curve and by the SMC-modified 
extinction curve computed by dust destruction numerical sim- 
ulations (SMC_40, Perna et al. 2003). We also indicate the posi- 
tion of the photometric filters for a system located at z — 0.695, 
the redshift of GRB 020405. Extinction curves that provide an 
acceptable fit to the data are emphasized with thicker lines. 



and normalization fluxes given by Fukugita et al. (1995). A 
7% systematic error was added quadratically to the magnitude 
errors to account for mag-to-flux conversion uncertainties and 
for intercalibration errors from different telescopes (see also 
Masetti et al. 2001). The resulting magnitudes are summarized 
in Tabled 

A simultaneous fit of the NIR-to-X-rays spectrum was per- 
formed with the X-ray spectrum corrected for the total pho- 
toelectric absorption and the optical SED corrected for the 
Galactic extinction but not corrected for an extra-Galactic com- 
ponent. For the intrinsic SED, we considered both a single pow- 
erlaw and a broken powerlaw, the latter with break energy be- 
low the optical range. In order to estimate the host galaxy ex- 
tinction we multiplied the continuum spectral model by an ab- 
sorption component. We tested different ISM dust composition 
and dust to gas ratios by assuming several type of extinction 
curves (see §3). 



3. Host galaxy extinction 

Since we do not know a priori the host galaxy extinction curve, 
we tested several types of environments (Fig.[Q with different 
dust grain size distributions, dust compositions and dust to gas 
mass ratios. The ratio between the hydrogen column density 
Nh and the visual extinction Ay provides a measure of the dust- 
to-gas ratio for a given dust distribution and composition. For 
simulated time-dependent extinction laws, the dust-to-gas ratio 
at the relevant time is derived from the simulations. 
We tested the following extinction curves: 

- the Galactic extinction curve (hereafter G) from Cardelli et 
al. (1989), for which AW^y = 0.18 x 10 22 ctrT 2 (Predehl & 
Schmidt 1995). 

- The Small Magellanic Cloud extinction curve (hereafter 
SMC) from Pei et al. (1992), for which N„/A v = 1.6 X 
10 22 cnr 2 (Weingartner & Draine 2000). 

- The attenuation curve derived for a sample of local star- 
burst galaxies (hereafter C) from Calzetti et al. (1994). No 
Nn/Ay has been measured for this case due to the complex- 
ity of the geometry of the dust and star distribution inside 
these galaxies (Calzetti et al. 2001). 

- Two extinction curves obtained by Maiolino et al. (2001) 
for the environment of a sample of AGNs. The Ql ex- 
tinction curve is computed assuming a power law dust 
grain size distribution dn(a) oc <r q da in the range a min 
0.005// < a < a max = 10// and q = -3.5. In this case 
N H /A V = 0.7 x 10 22 ctrT 2 . The Q2 extinction curve is 
derived assuming a ml „ = 0.005// < a < a max - l/i and 
q = -2.5 and it yields N H /A V = 0.3 x 10 22 cm" 2 . 

- Extinction curves resulting from numerical simulations 
performed with the code of Perna & Lazzati (2002). This 
code computes the temporal evolution of the dust grain 
size distribution, taking into account grain erosion due to 
UV and X-ray illumination. The ionizing continuum is tai- 
lored to the case of GRB 020405 (see above). Only the 
40 seconds of the burst emission were considered, since 
the overall fluence of the later afterglow was smaller than 
that of the prompt emission by at least an order of mag- 
nitude. Initial conditions were either a Galactic (final curve 
labelled G_40) or SMC (final curve labelled SMC.40) ISM. 
A uniform cloud with radius R = 7 x 10 22 cm and density 
n — 4 x 5 cm -3 was considered. We checked that changing 
this set-up would not change our conclusions. 

4. Results 

We fit the NIR/optical and the X-ray spectral energy distribu- 
tion, extrapolated to a common epoch (see §2), with a dust- 
absorbed power law model. The NIR-to-X-ray spectral index 
a was let free to vary. We found that the optical/NIR data 
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Fig. 2. NIR-to-X-ray absorbed power-law fits to the data of 
GRB 020405. A Q2 absorption curve is used in the upper panel, 
while a SMCL40 curve is used in the bottom panel. Optical data 
are shown both uncorrected (asterisks) and corrected (squares 
and triangles) for the host galaxy extinction. See Table[2]for the 
best fit parameters. 

corrected for the "standard" extinction curves, namely the G 
and the SMC, are incompatible with X-ray data (see Tab. |2j. 
Unacceptable fits were obtained also with the Calzetti, the Ql 
and the G_40 curves (see §3). 

Acceptable fits with comparable probability P(> x 2 ) of 
15%, were instead obtained with the Q2 and the SMC_40 
curves. These curves have the weakest dependence on wave- 
lengths among those considered (see Fig.[Q. We note that the 
relatively high^ 2 levels are due to the apparent bend of the con- 
tinuum in the NIR which cannot be reproduced by a power-law 
(Fig.E). 

We tested therefore a smoothed broken power law spec- 
trum (with the sharpness parameter s equal to 1, Granot & Sari 
2002) with break energy between the optical and NIR bands. 
The optical-to-X-rays spectral index was let free to vary within 
the 90% confidence level interval 0.8-1.0, found from the X- 
ray analysis. Assuming v e < vx, vo, from the X-ray spectral 
slope we derive an electron spectral index p = 2.0 + 0.4. We 
then compute the expected NIR spectral slope (corresponding 
to (p - l)/2 for v Nm < v e , Sari et al. 1998) that we let free to 
vary within the 90% confidence level interval (0.3 - 0.7). Even 
allowing for this extra degree of freedom, we found that the G, 
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Fig. 3. Same as Fig.|2]but with a smoothed broken power-law 
continuum model (see Table[3}. 

SMC, C, Ql and G_40 absorption models are inconsistent with 
the data (see Tab.[3j. 

For the models Q2 and SMC_40 we find an improvement 
in the fit with a 2.5<x statistical significance, according to the 
F-test (Fig. [3ji. They both fit the data successfully (Tab. [3}. 
We then checked the consistency of the measured Ay val- 
ues with the Nh z measured from the X-ray analysis and the 
Nn/Ay relationship expected (see §3). For the Q2 model, we 
found Nh = 0.7 x 10 22 cmT 2 assuming a power law model and 
Nh = 0.8 x 10 22 cirT 2 assuming a broken power law. These Nh 
values are consistent with the N Hz [(0.8 + 0.2) x 10 22 cm" 2 ] 
measured from X-ray analysis. For the SMC_40 model, even 
though the shape of the extinction curve fits the data well, a 
gray extinction with Ay > 2 is obtained only for Nh ~ 10 24 
cm" 2 , largely exceeding the X-ray measured Nh z - Formally a 
self consistent fit can be obtained assuming a dust-to-gas ratio 
100 times larger than the SMC one. This would imply a > 10 
times solar metallicity. We consider such physical conditions 
too extreme, even tough GRB (star formation) environments 
are expected to be especially dusty. 

5. Discussion 

We have studied the broad-band spectrum of the afterglow of 
GRB 020405 from NIR to X-ray bands. The optical-NIR and 
X-ray spectral indices are mutually consistent. However, this 
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burst is peculiar since it has an optical-NIR flux a factor of 
> 10 times dimmer than the extrapolation of the X-ray spec- 
trum. This can be attributed to an inverse Compton compo- 
nent in the X-ray band (Sari & Esin 2001). Alternatively, the 
optical brightness and spectrum could be affected by dust ex- 
tinction. We investigated the latter possibility. We found a self- 
consistent model that allowed us to constrain the properties of 
the dust distribution and its formation history. 

From a simultaneous NIR-to-X-ray spectral analysis, we 
found that the optical data can be fit with the external shock 
model only if corrected for an extinction curve weakly depen- 
dent on wavelength. Such an extinction curve can be the result 
of small dust grains coagulation and/or dust grain destruction 
(e.g. Maiolino et al. 2001, Perna et al. 2003). In Figure^ we 
plot the extinction curves derived from these two different pro- 
cesses (see §3). The remarkable similarity of these extinction 
curves, computed assuming different processes, is evident. 

The dust coagulation rate increases with density oc n l/2 
(Draine 1985). It is therefore favored in dense environments 
such as the cores of star formation regions where long GRBs 
are expected to explode. Incidentally, extinction curves weakly 
dependent on wavelengths have been inferred also for other 
extra-galactic objects, such as AGNs (Maiolino et al. 2001). 

A dust distribution biased towards large grains can also 
be due to grain destruction mechanisms (Waxman & Draine 
2000; Fruchter et al. 2001; Draine & Hao 2002; Perna & 
Lazzati 2002). We have attempted to reproduce the obser- 
vations starting from known dust distributions in the local 
Universe (Galactic and SMC). The evolution of the initial dis- 
tributions is followed as a result of the grain interaction with the 
burst ionizing flux (Perna et al. 2003). We find that the modified 
Galactic distribution cannot reproduce the observations, while 
the SMC one can reproduces them if an extremely large dust- 
to-gas ratio and metallicity are assumed. Since these conditions 
are hardly realized in nature, a scenario in which the large grain 
distribution is intrinsic to the host galaxy is favored. 

The difficulty in reproducing the results as a consequence 
of dust destruction can be understood since dust destruction re- 
duces the opacity besides modifying the reddening law. Since 
this burst requires a large extinction on a day timescale, it is 
likely that the opacity is provided by dust at a large distance 
from the burster, which is not affected by the burst photons. 
Prompt optical-NIR observations are better suited for the de- 
tection of dust destruction, since the evolution is caught "on the 
act". This would allow us to detect dust before it is destroyed 
and, based on the properties of what is left, infer the geomet- 
rical and physical conditions of the host galaxy ISM (Perna et 
al. 2003). 
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